EXAFS indication of double-well potential for oxygen vibration in Bai_ x K x Bi03 
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X-ray absorption spectra of oxide systems Bai-^K^BiOs and BaPb03 above Bi- and Pb-Z/3 absorp- 
tion edges were investigated. It was shown that oxygen ions move in double-well potential and their 
oscillations are correlated with charge carrier movement. Observed breathing-like oxygen vibration 
in double-well potential with large amplitude and low frequency causes the strong electron-phonon 
coupling and high T c values in doped BaBiC>3. Based on the experimental data, the model of rela- 
tionship of electronic and local crystal structures is proposed that is in a good agreement with the 
results of transport measurements, inelastic neutron and electron scattering, Raman scattering, and 
photoemission spectroscopy. In the framework of the model the possible reasons of superconductivity 
in perovskite-like oxides are discussed. 
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I. INTRODUCTION 

Though superconductivity in BaPbi-^Bi^Os (BPBO) 
was discovered significantly earlier |jj than in cuprates, 
the question of the nature of superconducting state in 
this oxide as well as in cognate system Bai-^K^BiOa 
(BKBO) is still unsolved. 

The structures of crystal lattice of copper oxide high 
temperature superconductors (HTSC's) and bismuth- 
based oxides have some important common character- 
istics. The both oxide classes have perovskite-like lattice 
with CuO„ (n=4, 5, 6) or Bi(Pb)C>6 complexes joined 
by the common oxygen ions. In bismuthates, the in- 
tersection of octahedral complexes in the three crystallo- 
graphic directions determines their three-dimensional cu- 
bic structure. The CuO„ complexes are joined in Cu02 
planes, which makes the two-dimensional structure of 
copper-oxides. 

Because of strong hybridisation of covalent Bi(Pb)6s, 
Cu3d - 02p a bonds, the above mentioned complexes 
are the most tightly bound items of the perovskite-like 
structure. Therefore such important peculiarities of per- 
ovskite structure as lattice instability in respect to soft 
tilting mode of CuO„ or BiOg complexes (see for review 
M ) and highly anisotropic thermal factors of oxygen ions 
vibration [ft, which point out the large amplitude of ro- 
tation oscillations, are inherent to the both classes of su- 
perconducting oxides and cause anharmonic vibrations of 
oxygen atoms that may be described by movement in a 
double- well potential . These structural instabilities 
of perovskite-like lattice can be related with the transi- 
tion to superconducting state f|||§|]. 

The layered structure of copper oxide compounds, 
presence of several non-equivalent copper positions, and 
a number of different Cu-0 bonds complicate the local 
structure analysis. At the same time, the simplicity of 
cubic three-dimensional structure of BPBO-BKBO sys- 
tems makes the interpretation of experimental data easier 
to a great extent. Relatively low temperatures of super- 
conducting transition T c ~ 13 K in BaPbo.75Bio.25O3 |l]] 




and T c ~ 30 K in Bao.6Ko.4Bi03 {|, the values of su- 
perconducting gap 2A(0)/fcT c = 3.6 ± 0.1 for BPBO § 
and 2A(0)/fcT c = 3.5 ± 0.5 for BKBO §], and a sizeable 
oxygen isotope effect Q allow one to rely on standard 
BCS-theory of superconductivity, not excluding a possi- 
ble realisation of other mechanisms. The simpler elec- 
tronic structure of s — p valence band of BPBO-BKBO 
systems in comparison with d— p band of cuprates favours 
the establishment of relationship of crystal and electronic 
structures in these compounds. 

However, even for these relatively simple systems there 
is no agreement so far on a number of crucial aspects: 
crystal structure symmetry and lattice dynamics flH 
electronic structure fl2| , bismuth valence state | [13| - 
The most part of unusual properties of BaBi03-family 
compounds mentioned in early review []l7f are still unex- 
plained. 

In addition, the average structural data contradicts the 
local ones. Integral crystallographic methods point to the 
simple cubic lattice in BKBO (x > 0.37) jO]. In contrast 
to this, the EXAFS-analysis of four nearest spheres of 
bismuth environment ]3(i| ] reveals a local tilting of octa- 
hedra by 4-5°, and Raman spectra evidence for a local 
lowering of symmetry from simple cubic |jl| . 

EXAFS-analysis of the nearest oxygen octahedral envi- 
ronment of Bi in metallic BKBO compound with x = 0.4 
previously was made in harmonic approach |32|3^] and 
pointed out at least one unresolved problem. Tempera- 
ture dependence of Debye- Waller factor a 2 (T) found by 
Heald et al. |32"[ can be described in Einstein approxima- 
tion, which is suitable for harmonic systems, only if one 
takes into account the temperature independent factor 
of 0.0025 A 2 . The attempt to explain this fact due to 
some static disorder with doping was denied since a 2 (T) 
of the next Bi-Ba and Bi-Bi shells show the absence of 
any disorder. Similar weakly varying u 2 (T) dependence 
was found in p3| . Thus, it was observed that amplitude 
of low temperature Bi-0 vibrations in metallic BKBO 
compound with x = 0.4 is too high (more than twice 
larger than in metallic BaPb03 compound) |^2],|4| ■ 
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Besides, our careful measurements [ p4j of temperature 
dependent EXAFS-spectra of undoped BaBi03 which 
were also treated in harmonic approach showed abso- 
lutely abnormal increasing of a 1 values of both Bi-0 
bonds with temperature decreasing from 90 K. We ex- 
plained these anomalies by the influence of anharmonic 
rotational vibrations of oxygen octahedra which become 
pronounced at low temperatures [ p9| , the model of con- 
nection of superconductivity with rotation mode asym- 
metry was proposed |35) . 

The above contradictions in description for local struc- 
ture lead us to necessity of EXAFS analysis of the nearest 
Bi-0 shell in anharmonic approach. 

In the present paper we report the results of temper- 
ature dependent EXAFS-investigation of BPBO-BKBO 
systems firstly treated in anharmonic approach based on 
the idea of such an analysis for apical oxygen atoms in 
YBa2Cu307_a proposed by J. Mustre de Leon et al. ||(| 
We have analysed EXAFS function x(k), using the new 
program "VIPER for Windows" |Q, by construction of 
the model potentials of atomic vibrations, subsequent 
calculation of the pair radial distribution function, and 
calculation of the model x(fe). This new approach for 
the EXAFS-analysis of Bi-based oxides gives us an op- 
portunity to investigate the character of oxygen atom 
vibrations and, in combination with our previous local 
electronic structure studies by XANES (x-ray absorption 
near edge structure) spectroscopy Q , to understand the 
nature of the structural phase transitions in these sys- 
tems and to explain practically all above contradictions. 

In Sec. [n] we describe the experimental details and 
present the procedure of data treatment. The general re- 
sults are given in Sec. [II. In Sec. IV we relate the local 



crystal structure with local electronic structure, which 
is the base of our processing of experimental data. In 
Sec. [v] we discuss the possible superconductive mech- 
anism and demonstrate how the oscillations of oxygen 
atoms in double-well potential contributes to supercon- 
ductivity in BKBO. 



II. EXPERIMENTAL AND DATA ANALYSIS 

In this work we have investigated ceramic samples of 
BaPb0 3 and BKBO with x = 0, 0.4, 0.5 synthesised as 
described in p9| ] . The materials were examined by x-ray 
powder diffraction for phase purity. The samples were 
controlled by transport and susceptibility measurements. 
BaBiOa showed semiconductor-like behaviour, the com- 
positions BKBO with x = 0.4 and 0.5 showed typical 
metallic p(T) dependence and superconducting proper- 
ties with T c ~ 30 K and T c ~ 16 K, correspondingly. 
BaPbOa samples were metallic but not superconducting 
at all temperatures. 

For the XAS measurements, a crushed fine powder was 
precipitated onto a micropore substrate. The thickness of 
samples was about two absorption lengths at the chosen 




FIG. 1. Experimental EXAFS x{k)k 2 (left) and its 
Fourier transform magnitude and imaginary part (right) for 
Bai-zK^BiOa (x = 0, 0.4, 0.5) and BaPb0 3 at 7K. 



absorption edge. 

The x-ray absorption spectra were collected at D-21 
line (XAS-13) of DCI (LURE,Orsay, France) synchrotron 
operated at energy 1.85 GeV and the average current 
~ 250 mA of positron beam at the L3 edges ofBi (13040.6 
eV) and Pb (13426 eV). Energy resolution of the double- 
crystal Si [311] monochromator (detuned to reject 50% 
of the incident signal in order to minimise harmonic con- 
tamination) with a 0.4 mm slit at 13 keV was about 2- 
3 eV. The low temperature measurements were carried 
out using a liquid helium circulation type cryostat with 
a temperature control of ±1 K. 

The background in the experimental spectra was re- 
moved as described in pQ] , taking care to remove the 
low frequency oscillations. EXAFS-functions x(fc)fc 2 ob- 
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tained from absorption spectra were Fourier transformed 
in the wave number range k from 1.5 to 16.5 A -1 , using 
Kaiser-Bessel windowing function. Back Fourier trans- 
form was done using a Hanning window from ~ 1 to 
~ 2 A corresponding to the first Bi-0 near-neighbour 
shell. In such a case, the number of independent ex- 
perimental points [[l2| was N exp — 2AkAr/ir + 2 w 11. 

The model EXAFS-function for pair atomic absorber- 
scatterer oscillations is constructed as follows. Suppose 
we know the potential of these oscillations as a para- 
metric function of interatomic distance. Solving the sta- 
tionary Srodinger equation for the particle with reduced 
mass of the atomic pair j37fl , one obtains a pair radial 
distribution function (PRDF) of atoms in i-th sphere: 



,(r) = A^|*„(r)| 2 e 



-E„/kT 



(1) 



where Ni is the coordination number, E n and ^ n are 
n-th energy level and its corresponding wave function. 
Given PRDF's, the model EXAFS function calculated as 



X(k) 



Mk)]/r 2 dr, (2) 



where k = ^j2m e /h 2 (E — E t h) is the photoelectron wave 
number referenced to the ionisation threshold E t h, r m i n 



and 



are determined by the windowing function of 



back Fourier transform. The phase shift 4>i(k) and scat- 
tering amplitude Fi(k) were calculated using FEFF-6 
code |Q for six-shell cluster with crystallographic data 
from neutron diffraction study [Tl| ] and using default set 
of the FEFF-6 parameters. The potential parameters 
were extracted from the model-to-experimental EXAFS- 
function fits. 



III. GENERAL RESULTS 

In Fig. [j] (left) we show the experimental x(k)k 2 for 
BKBO x=0, 0.4, 0.5 and for BaPb0 3 measured at L 3 
Bi(Pb) absorption edge at 7K. Good signal-to-noise ratio 
seen even for maximal wave number values k > 16 A -1 
indicates the high spectra quality. The absence of signal 
on the Fourier transform in the low-r range Fig. ^ (right) 
testifies for correct background removal procedure. 
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FIG. 2. Sketch of crystal structure of Bai-^K^BiOs in 
metallic (x > 0.37) and insulating (x < 0.37) states. 




FIG. 3. Experimental and model x(k)k 2 EXAFS for 
BaBiOa at 7 K for the first Bi-O near-neighbour shell (a) and 
the model potential with corresponding PRDF and energy 
levels (b). Uo is the energy of potential barrier 



A. BaBi0 3 

EXAFS-researches fl|,|§ of BaBi03 confirm the re- 
sults of crystallographic works |l^,|l(J[nj . According 
to them, there exist two inequivalent bismuth positions 
characterised by two Bi-0 bond lengths. The equality 
of coordination numbers of the two BiOg spheres points 
out that the BaBi03 structure represents the ordered al- 
ternation of small and large BiOg octahedra in barium 



lattice. Such an alternation together with static rota- 
tion distortion around [110] axis produce the monoclinic 
distortion of cubic lattice @|^|LlJ (Fig. |). As will be 
shown in Sec. to the larger soft octahedra corresponds 
the configuration BiOg, and to the smaller rigid octahe- 
dra corresponds BiL 2 Og. Here, L 2 denotes the hole pair 
in antibonding Bi6s02p CT » orbital of the octahedral com- 
plex. 
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In Fig. |(a) the experimental x(fc)fc 2 EXAFS for 
BaBi03 at 7K for the first Bi-0 near-neighbour shell is 
shown. The pronounced beating near 8 A -1 evidences 
for existence of at least two different lengths of Bi-0 
bonds. In all previous EXAFS-researches §l|,||,|2|-|3|] 
the EXAFS-function was quite successfully fitted in har- 
monic approximation as a sum of two independent har- 
monic functions with different bond lengths and Debye- 
Waller factors. However, temperature dependencies of 
Debye- Waller factors contradict to the harmonic Einstein 
model |34j . This argues against the independence of oxy- 
gen atom vibrations in two neighbouring octahcdra. 

In this work, we model the oscillatory Bi-0 potential 
as follows. Inequivalence of the two types of BiC>6 octa- 
hedra is due to presence or absence of a hole pair in the 
hybridised molecular orbitals Bi6s02p CT » (see Sec. IV). 
Suppose, the movement of oxygen atoms may transfer 
a hole pair from one octahedron to another. Such a 
movement exchanges the roles of two inequivalent octa- 
hedra and requires a double-well form of oscillatory Bi-0 
potential. Here, we take a parabolic form of each well 
Ui = Ki{r — ri) 2 /2 and U 2 = n 2 (r — r 2 ) 2 /2 which are 
joined continuously. Given the calculated x(k), defined 
by Eqs. ([!]) and (||), we perform a least-squares fit be- 
tween the model and experimental x(/c) over the range 
k = 2-16A -1 (see Fig. ||(a)). The six parameters deter- 
mined by the fit were E t h, ri, t 2 , Ki, n 2 , and N. The 
number of parameters here is the same as for fitting in 
the harmonic approximation. 

The analysis of parameters of double-well potential al- 
lows one to draw the following conclusions on the oscil- 
lations of oxygen atoms in BaBi03. (i) At low tempera- 
tures there exist two explicit peaks on PRDF, that is why 
x-ray and neutron diffraction detect the static distortions 
in BaBiOa. (ii) Even at low temperatures, the tunnelling 
probability between the two wells is non-zero. As it will 
be discussed in Sec. such a tunnelling is equivalent 
to the dynamic exchange BiL 2 Og <-> BiOg and explains 
observed activation conductivity in BaBi03. (iii) The 
temperature rise leads to gradual increase of probabil- 
ity of interwell tunnelling and to the structural transi- 
tion to cubic phase found at 750-800 K |l3|]. (iv) Non- 
cquidistance of energy levels implies that oscillations of 
breathing- or stretching types with several frequencies ex- 
ist: loq = Ei — Eq, loi = E 2 — Eq, uj 2 = £3 — Eq, etc. The 
latter are gradually manifested at high temperatures, (v) 
The positions of minima of double-well potential in our 
model are the average positions. They are spaced widely 
at maximal octahedra tilting and spaced closely at min- 
imal octahedra tilting. In the latter case the probability 
of interwell tunnelling is maximal. In this sense the tun- 
nelling frequency ujq is bounded by soft rotation mode 
frequency. 
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FIG. 4. Experimental and model x(k)k 2 EXAFS for 
BaPbOa at 7K for the first Pb-0 near- neighbour shell (a) 
and the model potential with corresponding PRDF and en- 
ergy levels (b). 



B. BaPb0 3 

The experimental x(k)k 2 EXAFS for BaPb0 3 for the 
first Pb-0 near-neighbour shell at all temperatures rep- 
resents a sinusoid with no beatings and phase breaks 
(Fig. |j) and fitted well, using a single parabolic poten- 
tial. This means that all Pb-0 bonds in BaPb03 are 
equivalent and the breathing- type distortion is absent. 



C. Bai-^K^BiOa 

The experimental x(fc)fc 2 EXAFS function for BKBO 
(x = 0.4) at 7K for the first Bi-0 near-neighbour shell 
is shown in Fig. |g(a), solid curve. Also, the model calcu- 
lated in harmonic approach (as in |32|]3g| ]) is shown as a 
dash-dot curve. It is seen that in the range k > 12 A -1 
the harmonic approach fails. This argues for the an- 
harmonic Bi-0 vibration behaviour similar to BaBi03 
case also in superconducting compositions BKBO. To a 
possible tendency in the cubic superconducting materials 
towards the same type of distortions as in BaBiOs was 
pointed earlier |32 
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FIG. 5. Experimental and model x(fe)fe 2 EXAFS for 
Bao.6Ko.4Bi03 at 7K for the first Bi-O near- neighbour shell 

(a) and the model potentials (single-parabolic and dou- 
ble-parabolic) with corresponding PRDF's and energy levels 

(b) . Total PRDF is shown by solid line. 



FIG. 7. Experimental and model xifyk 2 EXAFS for 
Bao.6Ko.4Bi03 at 7K for the first Bi-O near-neighbour shell 
(a) and the model potentials (single-parabolic and polynomial 
of degree 4) with corresponding PRDF's and energy levels (b). 
Total PRDF is shown by solid line. 




FIG. 6. The sketch of oxygen vibrations in Bai-^K^BiOs 
in the case of different (left) and equivalent (right) neighbour- 
ing octahedra with corresponding radial oscillatory potentials. 



The K doping of BaBiC>3 leads to partial replacement 
of the larger soft octahedra BiC>6 by the smaller rigid 
octahedra BiL 2 06 (see details in Sec. IV). This causes 
a decrease and disappearance of static breathing and 
tilting distortions, but keeps the different rigidities of 
Bi-O bonds. Hence, the movement of an oxygen atom 
depends on that to which neighbouring octahedra this 
atom belongs. If the neighbouring octahedra are differ- 
ent, the oxygen atom oscillates in double- well potential as 
in BaBi03. If the octahedra are equal, the oxygen atom 
oscillates in single-parabolic potential as in BaPbC>3 (see 
Fig. The statistical weights of these two cases de- 



pend on the potassium content x, and are (1 — x) and 
x, correspondingly. Here, the force constants of the two 
parabolas are assumed to be equal since their indepen- 
dent varying does not improve the fits, and the number 
of fitting parameters equals 7. The resulting model curve 
is presented in Fig. |](a), dotted. It is seen from Fig. ||(b) 
that the total PRDF is unsplit, so it is not surprising that 
crystallographic measurements reveal the cubic structure 
(e.g. at T = 10K a = 4.2742(1)). Meanwhile, the 
first momentum of the radial distribution function de- 
rived from our model is in a good agreement with the 
diffraction data |nj, which validates our calculated am- 
plitudes and phases. 

In this paragraph we address to the important ques- 
tion of statistical grounds for the choice among several 
possible models. Consider, first, the statistical chi-square 
function (labelled just like EXAFS-function, but this is 
the different value) 



X 



exp / , 



(3) 



where M is the number of data points in the fit, N exp 
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TABLE I. Bi(Pb)-0 parameters resulting from the fit to EXAFS data for BaPb0 3 and Bai-^K^BiOs with x = 0, 0.4, 0.5. 
r and I are positions (A) of potential minima and corresponding PRDF's first momenta, k is the force constant (eV/A 2 ). The 
index corresponds to parabolic potentials, the indices 1 and 2 correspond to the two parts of double-parabolic potentials. 
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was introduced above, Ej are the individual errors in 
the experimental data points. The latter were calcu- 
lated as average of yl/Io + l/It, lo an d h being the 
intensities measured in the transmission EXAFS experi- 
ment. The x 2 value must follow the \ 2 distribution law 
with degrees of freedom v = N exp — P, where P is the 
number of parameters varied during the fit. That is x 2 
must be less than the critical value X„ of the x 2 dis- 
tribution with v degrees of freedom and the confidence 
level c. For our model in Fig. ||(a) v 2 = 11 — 7 = 4 
and x\ — 5-3 < X°' 95 = 9.5, but for single-Gaussian 
model vi = 11 - 4 = 7 and X \ = 16-8 > X^ 95 = 14.1. 
Thus, our model meets the x 2 -test while the simple har- 
monic model does not. It is quite natural that having 
increased the number of parameters we got decreased \ 2 
value. But what the gain should be? The comparison be- 
tween the two models can be performed on the basis of 
F-tcst. If the difference x\ — X2 i s physically meaningful, 
not due to presence of the noise, that is if the simpler 
model cannot describe some features in principle, this 
difference must not follow the X 2 distribution law with 
V\ — v<x degrees of freedom. Provided that xi follows 
the x 2 distribution with v 2 degrees of freedom, the value 
F{vi - v 2l i> 2 ) = (X1/X2 ~ 1)^2/(^1 - v 2 ) must not fol- 
low Fisher's ^-distribution with v\ — v 2 and v 2 degrees 
of freedom. That is F must be greater than the critical 
value F£ 1 _ I/ U2 of the F distribution with v\ — v 2 and v 2 
degrees of freedom and the confidence level c. Compar- 
ing the two models in Fig. [j|(a), we find F(3,4) = 2.89, 
which is equal to F^ 8 f. Therefore, we can claim with 
84% probability that we propose a better model than 
the simplest single-Gaussian. 

The complete results for Bi(Pb) octahedral (N = 6) 
oxygen environment in BaPbOa and Bai_ x K x Bi03 with 
x = 0, 0.4, 0.5 at various temperatures are listed in Ta- 
ble H The uncertainties in the EXAFS distances and 
force constants are less than ±0.4% and ±50%, respec- 



tively. (These increments cause an increase of less than 
10% in the value of misfit.) It is seen from Figs. || and 
|5] that, in general, the positions of potential minima are 
not equal to the positions of PRDF's maxima. Because 
of that, among other parameters, we give the values of 
PRDF's centers of gravity, i.e. values that can be defined 
crystallographically. For BaBiOa, the centers of gravity 
were calculated for the two peaks of PRDF separately. 

The parameters of potentials were obtained through 
the fit and, of course, depend on the form of the model 
potentials. To elucidate this influence, we constructed 
the model EXAFS function, using the polynomial of de- 
gree four: U = -n{r - r ) 2 /2 + f (r - r ) 4 + k 2 /(16£), 
where the last term is introduced to zero the minimum 
values (Fig. |?j). The number of fitting parameters here is 
the same as for double-parabolic potential. It turned out, 
that for the both potential forms the mean frequencies of 
Bi-0 oscillations are practically identical for all tempera- 
tures. Because of this, in the present work we use double- 
parabolic potentials with parameters of clear meaning 
(potential minima positions and force constants) rather 
than a polynomial potential with abstract coefficients. 
Since the model EXAFS function weakly depends on the 
shape and value of interwell energy barrier Uq, we can 
not determine the Uq value from EXAFS measurements 
exactly. 

It worth to notice that in the case of BaBiOs the above 
model of oxygen vibrations has some limitation due to the 
existence of the static rotation (tilting) distortion ~ 11°. 
Thus, the distances observed from EXAFS spectra treat- 
ment are measured along Bi-O-Bi zigzag curve, but not 
along the Bi-Bi ([100]) direction as in BKBO. This may 
result in some difference between nl and k2 values (see 
Fig. |(b) and Table §) and some deviation of calculated 
frequency from real ones in BaBiOs. 
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IV. RELATIONSHIP BETWEEN THE LOCAL 
CRYSTAL AND LOCAL ELECTRONIC 
STRUCTURES 

A. BaBiO :i 

The co-existence in BaBiC>3 of two different types 
of octahedra with two different Bi-0 bond lengths and 
strengths reflects the different electronic structures of 
BiC>6 complexes. 

Octahedral complexes represent the most tightly 
bound items of the perovskite-like structures because of 
strong covalence of Bi6s-02p CT bonds. The valence band 
structure of BaBiCK is determined by overlapping of Bi6s 
and 02p orbitals PlJl^l, and, owing to strong BiBs-C^p^ 
hybridisation, the octahedra can be considered as quasi- 
molecular complexes |49| . Each complex has ten elec- 
tron levels consisting of a pair of bonding Bi6s-02p CT , six 
nonbonding 02^, and a pair of antibonding Bi6s02p CT » 
orbitals. A unit cell, which includes two octahedra, has 
38 valence electrons (10 from two bismuth ions, 4 from 
two barium ions, and 24 from six oxygen atoms). How- 
ever, the numbers of occupied states in the two octa- 
hedtral complexes are different: octahedron BiL 2 Og car- 
ries 18 electrons and has one free level or a hole pair L 2 
in the upper antibonding Bi6s02p CT » orbital, in octahe- 
dron BiOg with 20 electrons both antibonding orbitals 
are filled (Fig. ||). It is quite natural that BiL 2 06 oc- 
tahedra have stiff (quasi-molecular) Bi-0 bonds and the 
smaller radius, and BiOg octahedra represent non-stable 
molecules with filled antibonding orbitals and the larger 
radius. Because the sum of two nearest octahedra radii 
overcomes the lattice parameter a, the octahedral sys- 
tem must tilt around [110] axis, producing a monoclinic 
distortion in BaBi03 (see Fig. ||). 

The assumption of equal electron filling for nearest oc- 
tahedra (BiL 1 06+BiL 1 Oe) contradicts to experiments, 
since in this case equal Bi-0 bond lengths and local mag- 
netic ordering should be observed. 

Therefore our new scheme of bismuth disproportiona- 
tion 2BiL 1 06 — > BiL 2 06+Bi06 is in full agreement with 
charge balance, presence of two types of octahedron com- 
plexes and absence of any local magnetic moment . 

Because Ba 2+ in perovskite-type lattice is bound by 
pure ionic bond, the electron density is concentrated 
mainly in octahedra volume Vq J27J. Hence, the energy 
of the highest occupied level Ef is related to the octahe- 
dron radius R and the number of valence electrons N in 
a unit cell as 

E f ~ h 2 N 2 / 3 /m e Vo /S ~ h 2 N 2 / 3 /m e R 2 7 (4) 

where h is the Planck's constant, m e is the electron mass. 
This qualitative relation leaves out of account the devia- 
tion of Fermi surface from sphere |l7j]. Nevertheless, the 
value of Ef, which transforms to the Fermi level Ep at 
spatial overlapping of equal octahedra in the crystal, is 



20 electrons 1 8 electrons 




Bi0 6 BiL 2 6 
FIG. 8. The scheme of electronic structure formation for 
two different octahedral BiOe complexes. 



in close connection with the octahedron radius and with 
the number of valence electrons. 

According to expression (||) , the local electronic struc- 
ture is connected with the local crystal structure. In the 
left of Fig. ^(a), the case of hypothetical simple cubic 
structure of BaBi03 is shown. Real monoclinic structure 
arises from combined breathing (alternating octahedra 
with different radii R) and rotation distortions, which, 
according to ([I]), leads to the lowering of the energy Ef 
of the highest occupied Bi6s02p (T » orbital in BiO@ oc- 
tahedron in comparison with energy E h of unoccupied 
Bi6s02p CT » orbital in BiL 2 Og octahedron. Hence, the 
repetitive sequence of empty level Eh on the background 
of fully filled valence band represents the electronic struc- 
ture of the ground state of BaBi03 (Fig. ||(a), right). In 
such a system there are no free carriers and conductivity 
occurs at hopping of carrier pair from one octahedron to 
another with the activation energy E a = Eh — Ef. In this 
process the movement of hole pair in real space leads to 
change of large octahedra to small ones and vice versa 
or to a dynamic exchange BiL 2 06 <-> BiOe, in full ac- 
cordance with vibrations in double- well potential showed 
in Fig. and causes the hole-type conductivity. In this 
case, the activation energy E a means the pair localisation 
energy. 

This picture agrees with results of investigation of tem- 
perature dependencies of conductivity and Hall coeffi- 
cient: n(T) — n exp(— E' a j ksT) fijl, where the value of 
hole concentration no = 1.1 x 10 cm~ 3 coincides with 
the concentration of unit cells, and the activation energy 
E' a = 0.24 eV, if one takes into account that in the case of 
two-particle activation conductivity, the number of hole 
pairs is equal to the number of BiL 2 Og complexes (con- 
centration of which is half n ), and the activation energy 
E a = 2E' a = 0.48 eV. It should be mentioned that the 
possibility of the two-particle (bipolaron) conductivity as 
a most probable mechanism of conductivity in BaBi03 
was pointed to earlier [ ft7f . 

Pair splitting and hopping of a single electron from 
one octahedron to another cost energy and lead to elec- 
tronic structure reconstruction of the both octahedral 
complexes. Such a splitting can be achieved under optical 
excitation, which was observed experimentally as a pho- 
toconductivity peak at the photon energy hv — 1.9 eV 
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FIG. 9. The scheme of relationship between the local crystal and local electronic structures, (a) The ground state of BaBiC>3. 
Left — in cubic phase, right — in monoclinic phase. Eh is the energy of localised free level, Ef is the energy of the highest 
occupied Bi6s02po-» orbital. The black and white circles with arrows denote electrons and holes, correspondingly, with opposite 
spin orientations. The occupied states of Bi6s02p valence band are marked by grey. At the top, the local structure of octahedral 
complexes is shown. The [100] direction is equivalent to Bi-Bi direction, (b) The ground state of BaPbOs. The unoccupied 
band above the Fermi level Ef is shown, (c) The ground state of Bai-zK^BiOs, x < 0.37. The splitting of free level Eh 
at spatial overlapping of BiL 2 06 octahedra is sketched. E a is the activation energy E a = Eh — Ef. (d) The ground state of 
Bai-zK^BiOs, x > 0.37. The formation of unoccupied band above the Fermi level Ef at the percolation threshold reaching is 
sketched. 



p7[ . The optical gap E g is the energy difference be- 
tween the excited and ground states. In the excited 
state, BaBi03 has the local lattice of equivalent octa- 
hedra Bil/Oe and possesses the non-compensated spin. 
Admittedly, it must have anti-ferromagnetic ordering, 
as the ground state in undoped cuprates La2CuC>4 and 
Nd2CuC>4. Such optical excitation leads to the local dy- 
namical lattice deformation observed in Raman spectra 
as a breathing mode ^TOcm" 1 of giant amplitude under 
resonant coinciding of photon energy of Ar + laser with 
Eg Using lasers with other quantum energies 

destroys this resonant effect and leads to abrupt decrease 
of breathing mode amplitude |^0| . 

Thus, the scheme proposed accounts for the nature of 
the two energy gaps in BaBiC>3. The activation energy 
E a appears in transport measurements and is connected 
with coherent delocalisation of hole pairs owing to dy- 
namic exchange BiL 2 06 <-> BiC>6. The optical gap is 
radically differs from a traditional gap in semiconductors. 



It corresponds to energy difference between the excited 
and ground states. 

B. BaPbO :j 

In BaPbC>3, each octahedron also has ten molecu- 
lar orbitals, nine of which are filled and one, anti- 
bonding Pb6s02p CT » , is unoccupied as in octahedron 
BiL 2 06, since lead has one electron less than bismuth. 
Thus, all the octahedral complexes represent the bound 
molecules PbL 2 Og with equal radii. At spatial overlap- 
ping of PbL 2 06 in crystal BaPbC>3, wave functions of 
Pb6s02p CT « states become delocalised, unoccupied lev- 
els Eh of neighbouring octahedra split into a free carrier 
band, and, merging with the top of valence band consist- 
ing of initially filled Pb6s02p (T * orbitals, make the half 
filled conduction band (Fig. |^(b)). It is important to no- 
tice that unoccupied Bi(Pb)6s02p CT * orbitals in the case 
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of absence of spatial overlapping in BaBiOs behave as lo- 
calised hole pairs, but transform to the conduction band 
at the spatial overlapping in BaPbOs. Because of this, 
BaBiOa appears to be a semiconductor of p-type, and 
BaPbC>3 — a semimetal of n-type. Since the radius of 
PbL 2 C>6 complex is slightly greater than that of BiL 2 06 
complex, the Fermi level Ep in BaPbC>3, according to ex- 
pression (|J), lies lower in comparison with Ef in BaBiOs, 
which agrees with results of photoelectron spectroscopy 
& 

Thus, the lead valence state PbL 2 06 in BaPbOs is sim- 
ilar to the bismuth state in BiL 2 06 complexes in BaBiOs. 
According to Fig. [|, oxygen atoms in BaPbOs oscillate 
in the parabolic potential without any exchange between 
equal PbL 2 Og octahedra. Therefore, there is no charge 
pair transfer in semimetallic BaPb03 and the itinerant 
electrons behave as usual Fermi liquid. 



C. Bai zKrBiOa 

The doping of BaBiOs by lead or potassium leads 
to decrease of integral structural lattice distortions and 
causes essential changes in both local crystal and elec- 
tronic systems. Firstly, we discuss the changes in elec- 
tronic structure of BKBO as a simpler case, and then will 
extend the obtained conclusions to the BPBO system. 

Since K + ion has one valence electron instead of two 
of Ba 2+ ion, the substitution of each two Ba 2+ for two 
K + ions produces an additional unoccupied level or hole 
pair in a Bi6s02p (T * orbital and modifies the BiOg com- 
plex to the BiL 2 Og. As a result, the ratio of the num- 
bers of BiL 2 06 and BiO@ complexes changes from 1 : 1 
in BaBiOs to (1 + x) : (1 — x) and equals 7 : 3 in 
Bao.6Ko.4Bi03 and 3 : 1 in Bao.sKo.sBiOs. The spatial 
overlapping of BiL 2 06 complexes appear, which, taking 
into account their small radii and rigid bonds, contracts 
the lattice. This leads to decrease and disappearance 
(at x > 0.37) of both static rotation- and breathing-type 
distortions. The lattice is forced to contract, despite the 
ionic radii of Ba 2+ and K + are practically equal. 

Structural changes are accompanied by essential 
changes in physical properties of BKBO: at i k 0.37 
the phase transition insulator-metal occurs and the su- 
perconductivity arises that remains up to the dopant con- 
centration x = 0.5 corresponding to the solubility limit 
of potassium in BaBiOs. The type of the temperature 
dependence of conductivity changes from semiconduct- 
ing to metallic one, Hall coefficient changes its sign, and, 
in the normal state, BKBO compound with x > 0.37 
behaves as a metal with n-type conductivity p5[ . 

These changes are well described in the framework of 
the above scheme (Fig. ^ for BKBO). At the low dopant 
amount (x < 0.37), the contraction of the larger (soft) 
and the stretching of the smaller (rigid) octahedra bring, 
according to expression (^), Eh and Ef energies close 
together and the activation energy E a = Eh — Ef de- 



creases. As dopant concentration rises, the number of 
BiL 2 06 octahedra increases, and at x > 0.37, in the lat- 
tice arise the three-dimensional chains of spatially over- 
lapped BiL 2 06 octahedral complexes their unoccu- 
pied levels split and form the conductivity band, which 
is equivalent to the percolation threshold reaching that 
determines the insulator-metal phase transition. Here, 
the itinerant electrons show Fermi liquid behaviour with 
Fermi level Ep as in BaPb03. However, in contrast to 
the picture in BaPb03, in the BKBO structure there 
exist complexes BiOe through which itinerant electrons 
can not move because all the levels in these complexes 
are occupied. The movement of electrons occupying the 
upper Bi6s02p CT . orbital is possible only in pairing state 
at the dynamic exchange BiL 2 Oe <-> BiOe, since the un- 
pairing costs energy consumption. But, in contrast to 
the case of BaBi03, the static rotation- and breathing- 
type distortions of the octahedra disappear, their mean 
radii become approximately equal (and equal to half the 
lattice parameter), and the pair localisation energy ap- 
proaches zero: E a = E h - Ef — > (Fig. ^|(d)). In this 
case, the delocalised carrier pairs can freely move through 
the octahedral system that explains the appearance of 
superconductivity in BKBO (x > 0.37). This scheme of 
the local electronic structure is in full agreement with the 
picture of local oxygen vibration in double- well potential 
(see Figs.|,||). 

The unpairing of electrons, as in BaBiOs, is possible 
under optical excitation and is manifested by a pseudo- 
gap observed in reflectivity spectra even in the metallic 
phase of BPBO-BKBO systems. For instance, in BKBO 
with x = 0.4 the pseudo-gap is about 0.5 eV [^2| and 
about 0.6 eV in BPBO with x = 0.25. 

Our model suggests that in the metallic phase two car- 
rier types are present: itinerant electrons and pairs of 
initially (in BaBiOs) localised carriers. The co-existence 
of the two carrier types was confirmed by experiments 
on investigation of conductivity, Hall effect, and thermo- 
electric power 1 17 , as well as by photoemission spectra 
p3| , [5C| and Raman spectra investigations |lj||53|]. These 
results are in a good agreement with observed zero-bias 
conductance [ |54"| . The change in carriers behaviour from 
localised to itinerant at doping of BaBiOs by potassium 
was pointed to in x-ray absorption spectra analysis at 
the O K edge §f]. The analysis of EPR spectra ||§j56 
showed the presence of the localised carrier pairs, which 
also was confirmed by the observation of two-particle tun- 
nelling in the normal-state of BKBO |p4[]. 



D. BaPbi ^Bi^Os 

Practically the same changes in electronic structure 
arise at the doping of BaBiOs by lead. The elec- 
tronic structure of octahedral PbL 2 Og complex is entirely 
equivalent to that of BiL 2 Og complex. Unoccupied lev- 
els of overlapping octahedra in BaPbOs form the con- 







duction band, as at overlapping of BiL 2 Og in BKBO. 
At the doping of BaPb0 3 by bismuth (0 < x < 0.37), 
the case is similar to BKBO [x > 0.37). However, in 
BPBO the following combinations of neighbouring octa- 
hedra are possible: PbL 2 6 -PbL 2 6 , PbL 2 6 -BiL 2 6 , 
PbL 2 6 -Bi0 6 , BiL 2 6 -Bi0 6 . This leads to different lo- 
cal shifts of Eh and Ef energies, depending on different 
octahedra neighbouring pairs. At the further increase 
of bismuth content, the number of filled BiC>6 octahedra 
rises, the tilting and breathing distortions enlarge, the 
pair localisation energy raises, the spatial overlapping 
of unoccupied levels in PbL 2 Og octahedra disappears, 
which destroys the metallic type of conductivity and the 
system becomes a semiconductor with p-type conductiv- 
ity. 

Thus, superconductive properties of BPBO in our 
model are connected with dynamic exchange PbL 2 06 <-> 
BiOg. Unfortunately, it is practically impossible to ob- 
serve the double-well-potential oxygen vibrations for su- 
perconducting BPBO compositions (0 < x < 0.35) by 
EXAFS study because of overlapping of Pb L 3 and Bi L 3 
edges. Though an analysis of EXAFS spectra for these 
compounds is possible under elaborate treatment proce- 
dure p7| , the precise enough values can be obtained only 
for interatomic distances, but not for amplitude factors. 
For this reason, we do not present the EXAFS data of 
BPBO (0 < x < 1) in this paper. 



E. Photoemission spectra anomalies 

The analysis of photoemission spectra has revealed a 
serious contradiction concerning the experimental obser- 
vation of the splitting of the Bi4/(5/2, 7/2) spectral lines 
in superconducting compounds BKBO [£4j at the ab- 
sence of any peculiarities of those lines in parent BaBi03 
p3[ . The contradiction consists in that that though the 
analysis of x-ray diffraction and EXAFS data points to 
existence of two different Bi-0 bonds in BaBiOa, no split- 
ting of Bi4/ doublet lines was observed, which indicates 
the absence of considerable difference in bismuth valence 
states. The doping of BaBiOa by potassium relieves 
the monoclinic distortion and equalises the Bi-0 bond 
lengths up to a/2 . However, the Bi4/ spectral lines be- 
come broaden |23| ] and, at measurements on a high res- 
olution spectrometer, even split ]2^ j, which, in contrast 
with simple cubic lattice in BKBO (x > 0.37), points to 
existence of two different bismuth valence states. 

The above scheme of local electronic structure com- 
pletely resolves this issue. Indeed, in BaBi0 3 the bind- 
ing energies of the Bi4/ core levels in different BiL 2 06 
and BiOe octahedra are almost equal and hence, the 
Bi4/ doublet lines in photoemission spectra are unsplit 
(Fig. [lo](a)). Contrary to that, in Bao.6Ko.4Bi03, the 
binding energies of the Bi4/ core levels in the different 
octahedra E[, E' 2 and E\, E 2 differ on AE — Eh — Ep w 
E a , which causes the broadening (see curve x = 0.5 in 
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FIG. 10. Binding energies of Bi4/(5/2, 7/2) core lev- 
els in Bai-^K^BiOa: (a) for x = 0, right — experi- 
mental photoemission spectra from [^3| with energy resolu- 
tion AE — 0.45 eV; (b) for x — 0.4, right — experimen- 
tal photoemission spectra from Q| with energy resolution 
AE — 0.25 eV. Curve C is the result of convolution of the 
curve B with Gaussian function with 0.75 eV width. 




and splitting (see curve x = 0.4 
of the Bi4/ doublet lines. The 



Fig. |lfj(a) from 
in Fig. &Q(b) from 
value of the splitting coincides well with our estimation 
of E a = 0.48 eV, and the ratio of intensities / of the 
core-level peaks corresponding to the binding energies in 
different octahedral complexes BiOg and BiL 2 Og is in ac- 
cordance with the relative content of these complexes in 
Bao.6K .4Bi0 3 : I{E' 1 )/I(E 1 ) ~ I{E' 2 )/I{E 2 ) ~ 3/7. 



V. ON POSSIBLE NATURE OF 
SUPERCONDUCTING STATE 

Essentially anharmonic character of oxygen ion oscil- 
lations in soft rotation mode is the base of anharmonic 
model of the high temperature superconductivity 0. It 
was shown that anharmonic coupling constant \ s exceeds 
the constant in harmonic approach X p h'- 



As/Ap/j 



J*d 2 u/J 2 p h(u 2 )LU s » 1, 



(5) 



which explains high critical temperatures in cuprate 
HTSC compounds as well as in BPBO-BKBO systems. 
Here, d and lo s are the amplitude and frequency of oscil- 
lation in a double- well potential, {u 2 ) is the mean-square 
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displacement of the ions with the mean frequency ui in 
harmonic approach, J 2 h and J 2 are, respectively, the de- 
formation potentials of harmonic and anharmonic oscil- 
lations averaged on the Fermi surface. 

In that model, the oscillatory movement of oxygen ions 
in rotation mode with a large amplitude in the direction 
perpendicular to the Bi(Pb)-0 bonds was considered. A 
considerable gain in value of the coupling constant in ex- 
pression (||) is due to low frequency of soft rotation mode 
(lu s < lu) and the excess of its amplitude over amplitude 
of harmonic oscillations (d 2 ^> (u 2 )). 

Meanwhile, previous |lj,|27],|2^] as well as recent ( ]5§| ] 
and references therein) electronic structure calculations 
showed that such a movement hardly affects the Bi-0 
bond lengths, keeping sp(a) nearest-neighbour interac- 
tion nearly constant. Thus, despite the existence of 
double-well potential in rotation (tilting) mode, Mere- 
galli and Savrasov f^ ] obtained very small anharmonic 
contribution to A. Also, they assumed the breathing-type 
oxygen vibrations to be harmonic with small amplitude 
and high frequency and found too small A w 0.3 to ex- 
plain high T c in BKBO. 

Our investigation firstly showed that in BKBO the 
double-well vibrations with strong deformation poten- 
tial also exist along the Bi-O-Bi direction. Such vibra- 
tions have the breathing-like character and low frequency 
bounded by soft rotation mode. Thus, our results resolve 
the above issue and explain the reason for strong electron- 
phonon coupling in Bi-based oxides. Also, to these vi- 
brations we can assign the low- frequency part, below 
40 meV, of phonon density of states observed in neutron 
scattering for superconducting BKBO |39| and absent in 
calculations j58j. Besides, the anomalous phonon soften- 
ing along [100] direction |(3(J] becomes clear due to our 
observation of coherent breathing-like vibrations along 
[100] axis. 

Hardy and Flocken || evaluated values of A, using ab- 
stract model double-well potentials: 

FS 

A(T) = N(0)J2 £ ^Ze! 1 & - ^ (6) 

kk' n'>n 

where N(0) is the density of electron states at the Fermi 
level, M^k' is the e-ph matrix element between elec- 
tronic states \k) and |fc') on the Fermi surface, |n) and 
\n') are the oscillatory states with energies E n and E n i , 
/„ and /„/ are the thermal weighting factors; /„ = 
exp(-E n /fcT)/£ n , exp(-E n ,/kT). 

Given the oscillatory energies from our EXAFS- 
experiment, we have calculated A in terms of arbitrary 
multiplicative constant which is proportional to N(0) 
and Mkk 1 (Fig. [n]). Although "phonon" part of A is the 
strongest for BaBi03, this composition is not supercon- 
ducting. On the one hand, because of pair localisation 
energy E a . On the other, the rotation oscillations of rigid 
BiL 2 06 octahedra and ordinary phonons (of stretching 
and bending types) differ in collective character of mo- 
tion in Bi02 planes. In BaBiOa, rigid BiL 2 06 octahedra 
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FIG. 11. A's in terms of arbitrary multiplicative constant 
(solid lines) and mean oscillatory frequencies (dashed lines) 
for Bai-ajKajBiOa with x — 0, 0.4, 0.5 for double-well poten- 
tials obtained from EXAFS spectra treatment. 

are separated by soft BiOe octahedra and are not spa- 
tially overlapped. The partial replacement of the larger 
soft octahedra by the smaller rigid ones at the potas- 
sium doping leads to decrease and disappearance of lo- 
calisation energy E a (Fig. ^(c,d)) and to spatial coher- 
ence of rotation oscillations with length of several lattice 
parameters, depending on the doping level x. For su- 
perconducting compositions, T c means the temperature 
up to which A is considerably decreased (Fig. [ll]) and/or 
the coherence in oscillations of neighbouring octahedra is 
thermally destroyed by ordinary phonons. 

The ground state of BaBiOa can be considered as a 
bipolaronic state as well, and the movement of carrier 
pairs correlated with oxygen atoms oscillations at the 
dynamic exchange BiL 2 06 <-» BiOe evidences for possible 
application of bipolaron theory pl[ ] , which is supported 
by the small size of the pair (octahedron size) and by the 
existence of the pair state above T c . 

VI. CONCLUSION 

From the EXAFS investigation of BPBO-BKBO sys- 
tems, the local crystal structure peculiarities connected 
with non-equivalence of BiL 2 06-Bi06 octahedra were ob- 
served. It was found that oxygen vibrations are well de- 
scribed using double- well potentials, which leads to the 
strong electron-phonon coupling due to coherent mod- 
ulation of Bi-0 bond lengths with low frequency and 
causes correlated carrier pair movement with oxygen ions 
oscillations. The underlying relationship between the 
local crystal and local electronic structures was estab- 
lished that explains the full list of unusual properties of 
BPBO-BKBO systems: the existence of two energy gaps 
(transport E a and optical E g ) in BaBiOa; the mech- 
anism of two-particle conductivity in BaBiOa; the co- 



ll 



existence of two different carrier types; pseudo-gap ob- 
servation in metallic compositions of BKBO and BPBO; 
the observation of localised pairs from EPR spectra; the 
observation of low-frequency (< 40meV) phonons; the 
nature of concentration and temperature phase tran- 
sitions; the contradictions between local (Raman, and 
EXAFS) and integral structural methods for metallic 
BKBO; the anomalies of XPS spectra in BaBiC>3 and 
BKBO. We observed the correlated movement of carrier 
pairs and low-frequency breathing-like oxygen octahedra 
vibrations, which corresponds to superconducting state. 
The model proposed combines some principal features of 
the real-space pairing Q , anharmonic models j|,|5| and 
bipolaron theory |n) of high-T c superconductivity. 

The likeness of rotation mode peculiarities of Bi06 
octahedra and CuO„ complexes and the anomalies in 
temperature dependencies of Debye-Waller factors in 
cuprates J62[ allow one to hope that similar model ap- 
proach can be applied to the Cu-based superconductors. 
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